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' Abstract 

This paper develops a modification of the dressing method based on inhomogeneous 
linear integral equation with integral operator having nonempty kernel. Method allows 
one to construct the systems of multidimensional Partial Differential Equations (PDEs) 
having differential polynomial structure in any dimension n. Associated solution space is 
not full, although it is parametrized by certain number of arbitrary functions of (n — 1) 
variables. We consider 4-dimensional generalization of the classical (2+l)-dimensional 
■ S'-integrable iV-wave equation as an example 

(N 

1 Introduction 

*Ch ■ 

Completely integrable nonlinear Partial Differential Equations (PDEs) became an attractive 
field of research after the discovery of complete integrability of the Korteweg-de Vries equation 
[1]. S'-integrable [2, 3, 4, 5, 6, 7, 8] and C-integrable [9, 10, 11, 12, 13, 14] systems are mostly 
remarkable among multidimensional nonlinear PDEs. The dressing method is one of the pro- 
moted methods for constructing and solving S'-integrable PDEs. After the original version of 
the dressing method [2] several important modifications have been developed [3, 4, 5, 15, 16, 17] 
(see also [7, 8]). All of them are based either on a Riemann-Hilbert or on a <9-problem [18, 19], 
which are uniquely solvable linear integral equations for some matrix function U (A; x) depend- 
ing on spectral parameter A, where independent variables of nonlinear PDEs appear as a set of 

additional parameters X — (^1 , X2 ,...). 

This paper is based on the modification of the <9-problem introduced in [20, 21]. It was 
shown [21], that classical S-integrable PDEs, some types of C-integrable PDEs and some types 
of "mixed" PDEs (i.e. equations having S- and C-integrable systems as the particular cases) can 
be studied by the dressing method based on the following integral equation, which is equivalent 
to the <9-problem [4, 5]: 



$(A;z) = J ^(X,fj,;x)U(fi;x)dQ(fx) =^(X,fx;x)*U(fj,;x) = $U. (1) 
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Here A and fi are spectral variables, U is the unknown matrix function. The matrix functions 
$ and \1/ are defined by some extra conditions, which fix their dependence on the additional 
parameters Xi. Q is some largely arbitrary scalar measure in the /x-space. Apart from Q, all 
the functions appearing in this paper are Q x Q matrix functions. 

To describe the classical completely integrable equations, kernel of the integral operator ^ 
must be trivial, i.e. dimker ^ = 0. It was shown in [22] that this requirement is too restrictive. 
More general dressing algorithm must be based on the integral operator with dimker\l/ > 0. 
So, some examples of nonlinear PDEs corresponding to dimker\l/ = 1 have been derived in 
[22]. Generalization dimker\l/ > removes any formal restriction on dimensionality of the 
constructed nonlinear PDEs, although derived systems have restricted solution spaces. Namely, 
n-dimensional nonlinear system admit at most (n — 2)-dimensional solution space, n > 1. 
The increase of dimensionality of nonlinear PDEs in comparison with classical S'-integrable 
equations happens due to the presence of the external dressing function in the algorithm, whose 
dimensionality has no restriction. Remember that the classical dressing algorithm uses only 
the internal dressing functions whose dimensionalities are strongly restrictive and usually equal 
1. However, the new version of the dressing algorithm has an important common feature 
with the classical algorithm. Namely, derivatives *ff Xj (\, x) are separated functions of the 
spectral parameters A and /i. This fact fixes dimensionalities of the internal dressing functions 
mentioned above. This property has been removed in [23] (i.e. ^ Xj (X, /j,; x) are not separated 
function of spectral parameters there) keeping dimker^ = 0. But there is serious obstacles for 
construction of the explicit solutions, which have not been found in [23]. 

Another problem appeared in [22] is a complicated "block" structure of the derived nonlinear 
PDEs, which is not typical for the equations of mathematical physics. 

Combination of basic ideas of both refs. [22] and [23] enriches the solution space and im- 
proves structure of the system of nonlinear PDEs derived in this paper. The dimensionality n of 
the resulting system of nonlinear PDEs has no restriction and its solution space is parametrized 
by certain number of arbitrary functions of (n — 1) variables, although the full integrability is 
not achieved yet. Nevertheless, this is a progress in comparison with the results of [22], where 
the solution space of the derived system of nonlinear PDEs is parametrized by arbitrary func- 
tions of (n — 2) variables. Structure of nonlinear PDEs is also improved: system of equations 
has differential polynomial structure. 

Algorithm represented in this paper is based on the following modification of the eq.(l): 

Q . 

$( sfc )(A;x) = / T {n) y {s) (Kv,x)U {nk) (n;x)dn(n) = (2) 

n=l 
Q 

J2T^¥ s \X,fi-x)*U {nk \fi;x) = 

n=l 

Q Q 

y (sn) (K A«; x) * U {nk) (fi; x) = J2 ¥ sn) U {nk \ 

n=l n=l 

^(A, fi; x) = T^¥ S \X, n;x), s = l,...,Q, 

which must be solved for the matrix spectral functions u( nm \ Here are constant diagonal 
matrices, functions $^ and \&( s ) are called the internal dressing functions and will be specified 
below. We write superscripts inside of parenthesis in order to distinguish them from the power 
notations. 
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It was observed that the integral equation (2) may be simplified taking 



Ti n) = W = { i' n T (3) 

and the diagonal function $0 m ) in all examples studied in this paper. So that the integral 
equation reads: 

< } <MA; x)=J2 e^S(A, ^ x) * U<f =► (4) 

n,7=l 

" E?=i^(A^;x)*^f (^;x) = < ) (A;x), a = 
£? =1 *S(A, s) * x) = 0, a ? /? ' U 

Hereafter (except Sec. 3) double subscript (usually Greek) means the element of the appropriate 
matrix and single subscript denotes the nonzero element of the appropriate diagonal matrix 
unless different is specified. 

At first glans, following the philosophy of the dressing method based on the integral equa- 
tion, integral equations for (A; x) and U^ k \X;x), a ^ (3 are decoupled. This might result 
in appropriate decoupling of the system of nonlinear PDEs. Moreover, the integral eqs. (5a) 
and (5b) with different values of the indexes a and (3 might be decoupled as well. However, 
we will show that general situation is different and the above statements are partially correct. 
Namely, we will obtain an example of the complete system of nonlinear PDEs written for the 
fields produced by U^\ k \X;x) and u!^q\X;x), fc,7 = 1,...,Q. However, more general PDEs 

may couple fields produced by all functions uj^ k \X; x), k, a, (3, 7 = 1, . . . , Q (see Sec. 2. 2, item 
(d.3) of Proposition). 

In the next section (Sec. 2) we represent a generalization of the dressing method for the 
classical (2+l)-dimensional A^-wave hierarchy. In Sec. 2.1 we discuss the case dimker\l/ = 
giving the classical (2+l)-dimensional A^-wave equation 

BV } - BWvtB® + [v yi ,B^] - [v y2 ,B^] + [[v, flW], [v, B^]] = 0, (6) 

where B^ are constant diagonal matrices, v is Q x Q matrix field, t is time and y^ i — 1,2, 
are space variables. Its generalization, system of 4-dimensional 1-st order quasilinear PDEs 
solvable by our algorithm with dimker^ = 1, is proposed in Sec. 2. 2. The whole system of 
nonlinear PDEs derived there may be separated into two subsystems. The first subsystem is 
the complete system of evolution equations for the Q x Q matrix fields p = 1,2, v 
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and u: 

*Sr ;4) ^S + E C^^S - E ^^isr 5 + ( 7 ) 



m=l 7=1 



Q Q 2 



E $ vCr + E E ^>S^ 4 " i0) = °> 



7,4=1 7=1 j 0= l 

4^/3 7 ^a 



# 9i4) a«$ + E *i 4;9;m) ^$ - E ^M? wu) + (8) 

m=l 7=1 

E $Se T «T ] + E E &MS )T # m * ) = °> 



7,5=1 7=1 i =l 

7#« 



S^dtVaf, + S ap m)d ymVaP ~ Yl V *l V lP T a-7 V) + ( 9 ) 
m=l 7=1 

7?W/3 

Q Q 2 

E u «i v w T «S v) + E w ^ E ^7/3 )T «7r" ;4o) = ' a ^ & 



7,4=1 7=1 in=l 



af^fltti^ + S^ 1 ^^^ - £ V^flQ™* + (10) 



7 = 1 



Q 2 



E «* y ^Sr ) + E E ^ T £ m) = °> 



7,4 = 1 7=1 ig = l 

4^a 

The equations of the second subsystem must be viewed as a symmetry (i.e. compatible) con- 
straints to the system (7-10): 

s^' 3) d y ^ + £ sf^dyM - E + (11) 

m=l 7=1 

E ^SS^St* + EE ^Tg™** = o, 

7,5=1 7=1 i =l 

75^a 
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s ( :p 3) d y3 v a , + £ s^d ym v aP - J2 + (12) 

m=l 7=1 

E wv^+E^E^Sr*^ . 

7. 4=1 7=1 «0 = 1 

Q 

s<*™%uafi + s^ l) d yi u aP - E v ai u^T^ + (13) 

7 = 1 
7^" 

E ^T^^r u) +EE ^^^if u ^ = o, j = 2, 3. 

Although eqs.(7-10) represent the complete system for unknown fields, the dressing method 
provides only such solutions which satisfy the symmetry constraints (11-13). Here the constant 
parameters T and s are expressed in terms of the dressing data, see item (c) of Proposition, 
eqs.(52) for details. Not all of them may be arbitrary. This system has the remarkable reduction: 
if u — then the eq.(9) yields the iV-wave equation (6) while the eq.(8) yields the following 
linearizable system 

drJll + E E $h { $ T tt qqq '' i0) = °> 3 = 3, 4, (14) 

7 = 1 i = l 

where r a are characteristic variables, see eq.(53) for more details. Linearization of this system 
will be demonstrated in Sec. 3.1. 

Note, that another system of PDEs having S- and C-integrable systems as particular reduc- 
tions has been derived in [21]. Similar to the PDEs derived in [21], the system (7-10) consist 
of 4 types of nonlinear equations, which defer by dimensionalities of their linear parts. It is 2 
in the eq.(10), 3 in the eqs.(8,9) and 4 in the eq.(7). However, the structure of the nonlinear 
PDEs is significantly different. In particular, some coefficients of the system (7-10) depend on 
the parameters R!yi \ 7, 5 = 1, . . . , Q, reflecting the fact that dimker^ = 1 (see eq.(41)). Such 
parameters do not appear in [21] since dimker\l/ = therein. 

Let us clarify which kind of initial-boundary value problem may be solved, in principle, 
numerically. It may be easily seen that the fields may be given arbitrary initial conditions 
(i.e. values at t — 0) in the whole three dimensional space (yi, yi, ])•&)■ The fields and v 
may be given arbitrary initial conditions only on the plane (for instance, 1/3 = 0) because of the 
constraints (11,12). In order to define their initial conditions in the whole space (yi, 1/2,1/3) we 
must solve the system (11,12). Field u is even more restrictive because of the constraint (13). 
Initial condition for u may be arbitrary on the line ( for instance, y% — y2 — 0). To define initial 
condition in the whole space (2/1, 2/2,2/3) we must solve the system (13). Finally, t-evolution of 
all fields may be established solving the system (7-10). We conclude that correctly formulated 
initial-boundary value problem requires the following initial-boundary conditions: two arbitrary 
matrix functions of 3 variables w^\ t =o, p = 1,2, three arbitrary matrix functions of 2 variables 
q^\t= y3 =o, p — 1, 2 and v\ t=y3= o, and single arbitrary matrix function of one variable u\ t= y 3=y2= o. 
If such initial-boundary data are provided by the dressing method we would consider the system 
(7-10) with compatible constraints (11-13) as a completely integrable one. 
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However, we will show that the represented dressing algorithm may supply only single 
arbitrary matrix function of three variables, two arbitrary matrix functions of two variables and 
single arbitrary matrix function of one variable. For this reason, using our dressing algorithm, 
we are not able to prescribe arbitrary initial data to all fields. This fact causes us to refer to 
such equations as partially integrable PDEs. It will be shown by construction that they admit 
infinitely many commuting flows. 

Remark, that another system of nonlinear evolution PDEs with symmetry constraints for 
some fields in the context of the dressing method has been derived in [24]. 

Detailed derivation of the system (7-13) is given in Appendix, Sec. 5. Sec. 3 describes the 
solution space to the 4-dimensional system (7-13). Conclusions are represented in Sec. 4. 

2 First order nonlinear PDEs: higher dimensional gen- 
eralization of the /S-integrable N-wave equation 

Dressing method for the (2+l)-dimensional iV-wave equation can be generalized to higher 
dimensions. Here we consider a variant of such generalization. For this purpose, function 
^^(A, fi; x) in (4) must be related with the dressing functions & sm \\] x) and C^ m \fi; x) by the 
following bi-linear system of equations introducing dependence on the variables Xj, i — 1, 2, . . . : 

Q 

d Xj ¥°\\,Kx)-A®d Xl ¥ a \\,r,x) = J2^ sk) ^^)B m C (k \^,x), (15) 



k=i 



B m = I, j > 2. 



Here / is the identity matrix, and are constant diagonal matrices, dressing functions 
C^(A; x), k — 1, . . . , Q, will be specified below. As far as the system (15) is overdetermined 
system of PDEs for the functions ^( s \ this system must be compatible. Compatibility condition 
reads 



[d Xj ($ (sk) (\;x)B^C (k \fi;x)^ - d Xi (V sfc )(A; x)B^C^(n; a;)) - (16) 

fc=i 

A U) dxi (A; x)B^C^\fi; a;)) + A^d Xl (& sk \\; x)B^C {k \^ x)J 



Eq.(16) consists of terms represented by products of functions depending on single spectral 
parameter. We would like to separate PDEs involving different spectral parameters, i.e. ei- 
ther parameter A or //. This is possible due to the diagonal form of $( STra ) which provides 
commutativity of A^ and $( sfc ); 

A«$( sfe )(A; x) = & sk \\; x)A®, V i, s, k, (17) 

We put i — 2 in (16) without loss of generality and separate PDEs involving different spectral 
parameters: 

d x ^ {sk) (X;x) = d Xl <S> {sk \\;x)P m + d X2 <5> {sk) (\; x)B {kj) , (18) 
d X] C^ k \^x) = P^d Xl C^(fi-x) + B^d X2 C^^;x) (19) 

P (kj) = A u) _ A m B m^ j>2 (20) 
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Eqs.(17) and (18) must be used for the derivation of nonlinear PDEs as follows. First, 
substitute ¥ sm ^ from (4) into (17): 

Q Q 

A^ sn \\, //; x) * U {nk \fi; x)-J2 ^ (sn) (A, fi; x) * U {nk \fi; x)A® = 0, (21) 

n=l n=l 

Second, substitute $ (sfc) from (4) into (18) (we use (17) to result in (22)): 
Q 

\ d Xj (^ sn \X, fi; x) * U^di; x)) - A^d xi (* (sn) (A, /i; x) * U^ nk \/j,; x)) - (22) 



n=l 



d X2 (¥ sn \X, n- x) * U^\n; x)) - A®d xi (¥ sn \X, x) * U (nk \^ x)))b™ 



0. 



Eqs. (21,22) must be reduced to homogeneous equations (see eq.(23)). For this purpose we 
substitute ^ X J from (15) into (21,22) which results in the next homogeneous equations: 

Q 

J2¥ sn \X,n;x)*L (ji ' nk \^x) =0, 2 = 1,2, (23) 

71=1 

where 

L^ nk \\-x) = U {nk) B (j ' n) , j>2, (24) 
L {j2 ' nk) (X; x) = d Xj U {nk) + d Xl U {nk) B^ kj ' n) - d X2 U^ nk) B {kj) + (25) 
Q 

[/(nil) ^ Jj(hj) v {nhk) _ ^(niifc) g(kj)^j ^ j > 3 ; 

ii=l 

fields v( hk ) are introduced by the formulae 

v^ lik \x) = C*(/i; x) * U {lk \^ x), (26) 

and the diagonal matrices B&n) and B^ kj ' n ^ are given by the formulae 

B U;n) = A^I-A ij) , B {kj ' n) = B^A® - A^I. (27) 

In addition, applying operators (d Xj — AW Xl ) to the eq.(23) with % — 1 and j = jo, one gets 
one more homogeneous equation: 

Q 

^ {sn) (A, fi; x) * L^ nk \n- x) = 0, (28) 

n=l 

where 

Q 

L U3 > nk) (\; x) = (d Xj U {nk) - d Xl U {nk) A^ + ^ jj {nil) B {hj) v {nhk) ^jB {kr ' n \ j > 2. (29) 

h=i 

Comparing the eq.(20) with the eq.(27) we see that 

p(kj) _ _jg{kj;a) ^ j-gg-j 

Although this paper is devoted to the integral equation (4) with dimker ty( sn ^ > 0, we 
consider the case dimker ^( sn ) = in the next subsection for the sake of completeness. The 
case dimker i]A sn ) = 1 will be considered in Sec. 2. 2. 
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2.1 dimker ^j/( sn ) = 0. ^-mtegrable (2+l)-dimensional iV-wave equa- 
tion 

If dimker ty( sn ^ = then the homogeneous equation 

Q 

^ (sn) H (nfc) = (31) 

n=l 

has only the trivial solution H*™^ = 0, i.e. 

L (n-M) ;= u { ^ k) B { i' n) = (32) 
:= d X] U^ ) + d Xl U^ ) B^ n) -d X2 U^ ) B^ ) + (33) 



i,7=l 



:= U^-^l^^ ( 34 ) 

V i,7=l / 

Since £% ;n = 0, eq.(32) tells us that = if n ^ (3, and consequently t^™^ = if n ^ (3. 

Then, eq.(34) is identical to zero, while eqs.(33) with different values of (3 become decoupled: 

V tj (X; x) - V t2 (X; x)B® - V(X; x) [v(x),B&] = 0, j = 3, 4, . . . , (35) 

where, for fixed (3, 

V ak = U% k \ v* = vffi\Bjp = B<M, a,i,fc=l,...,Q, (36) 
dt.=d x .-Afd xl , j = 2,4,.... 

Eqs. (35) represent the linear overdetermined system for the spectral function V. Compatibility 
condition of (35) yields the classical .S-integrable (2+l)-dimensional iV-wave equation: 

K,BW] - [v t .,BW] + B^v t2 B^ - BWvtB® + [[v, B^% [v, B^}} =0, k ^ j = 3,4, . . . .(37) 

This example justifies the fact that our dressing algorithm with dimker ^( sn ) = gives rise to 
the classical S'-integrable models. 

2.2 dimker ^( sn ) = 1. Higher dimensional nonlinear PDEs 

Let dimker \|>( sn ) = 1. Then the solution space of the homogeneous equation (31) is parametrized 
by the arbitrary Q x Q matrix functions f^ %k \x) ) i, k — 1, . . . , Q: 

Q 

(U h ) (nk \X; x) = J2 H {ni \X- x)f i - ik \x). (38) 
i=i 

This equation establishes a linear relation between any two solutions U^ nk \ L^ 1;nk \ L^ 3 ' nk \ 
j > 2 and L^ 2;nk \ j > 3 of the homogeneous integral equation (31). These linear relations 
represent a new overdetermined system of linear equations for the spectral functions 
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Details of derivation of the system of linear PDEs for the spectral function £/(™ fc ) as well as 
derivation of the assotiated system of nonlinear PDEs are given in Appendix, Sec. 5. Here we 
collect all basic results in the following Proposition. 

Proposition. Let dimker ^( sn ) = 1 in the integral equation (4) where the dressing func- 
tions ^ (s) , $ (sfc) and C (fc) are solutions of the eqs.(15, 17,18, 19). Then 
a) Matrix functions of x (in other words, fields) may be introduced by the formulae 

v (nik) ^ = C (i) ( A . x) + jj(nk) ( X . x ^ v (nik;l) ^ = ^ ^. ^ + jj(nk) ^. % ^ ^) 

w ( nk ) ( x ) = G(X;x)* U {nk) (A; x) , w (nk ' p) (x) — G Xp (X;x) * U (nk) (A; x), 



w 



(nk\ps) 



(x) = G XpXs (\;x)*U^ k \X;x),..., 



where n,i,k = 1,...,Q, p,s = 1,2,..., function G(X;x) is the external dressing function, 
whose role will be clarified below, see Sees. 3 and 5. These fields are related with each other by 
the system of nonlinear PDEs, see eqs. (101-103,106-111) in Appendix. If G is arbitrary, then 
the system of nonlinear PDEs may not be completed. 

b) Presence of arbitrary matrix functions f^ k \x), i,k = 1, . . . , Q in the solution space of the 
homogeneous eq.(31) (which is indicated in the eq. (38)) allows one to impose a largely arbitrary 
relation among the matrix fields (39): 



T nk (all fields ) = 0, n,k = l,...,Q, (40) 

where T nk are Q x Q matrices. The only requirement to this relation is its resolvability with 
respect to f^ k \x). 

c) Let relation (40) be taken in the next simple form 

w (nk) = R (nk) ( 41 ) 

(where R^ nk \ n,k = 1, . . . , Q, are constant non-degenerate Q x Q matrices), and the external 
dressing function G(X; x) satisfies the following system of linear PDEs 

G Xj (X;x) = V^G Xl (X;x)+V^G X2 (X;x), j > 2 (42) 

with arbitrary diagonal matrices J>^ 1 \ i — 1,2. Replace Xi by the new independent variables 

d tj = d Xj - A?d xi , j>2 => t = t 5 , Vj = t j+1 , j = 1, 2, 3. (43) 

Then the complete system of nonlinear PDEs can be derived for the following fields 



(p) _ (l/3;p) 



a/3 — w al 

7,5=1 



$ = E < 5;P) £? ;1) (^C k, a = 1, . . . , Q, p = 1, 2. (44) 



Vap^vM, Ua^^v^B^iR- 1 )^, i,p = l,...,Q. 

7,5=1 

This system of PDEs is naturally separated into two subsystems. The first subsystem is the 
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complete system of evolution PDEs: 



3 



E (56;l/3p) ._ (4;«,;4)o (p) , (4;«,;m)o Xp) _ \^ Jp) T (4;wwv) 



m=\ 7=1 

Q Q 2 

g Q7 "8{)-L al 6l3 + ^"7 "V J cry/3 — U ' 

7,<5 =1 7=1 io = l 

5^/3 7^" 

£(57;W ;= s (W) dtq W + £ _ £ + (46) 

m=l 7=1 

E AW-^ + ee ^ T $ qqq '' l0) = 

7,5=1 7=1 i =i 

75^« 



Q 



m=l 7=1 

7t>W/3 

Q Q 2 

7><5 =1 7=1 *0 = 1 

Et' la " ] ■= s^d t u a(3 + s^d yi u a , - J2 v a ,u^T^u) + (48) 



7=1 



Q 2 



E ^u 5P T^r u) + E E ^t^^ = 0' 

7,*=1 7=1 ?0 = 1 

The second subsystem must be viewed as a system of symmetry (i.e. compatible) constraints 
to the system (45-48): 

Emm := sf^d mq § + E sf^d y J:} - E ^u^ T S qwu) + ( 4 9) 

m=l 7=1 



E «S^Sr ) + E E ^tSt** = °> 



7<5=1 7=i i =l 
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E 



(42; lap) 
11 



((j+l)3;lc*/3) 



Q 



(3;«;3) o s ( 3 ^ m ),g . _ ?; 7; „T {3;wt,) + 

6 a/3 °y3 V otf3 -+- 6 a/3 U y m U al3 2^1 a ^ ^ "7/? + 

m=l 7=1 

77^/3 

E w^^+E^E^Sr*^ . 

7,<S=1 7=1 io = l 

Q 

7=1 

E w*r£r° +ee ^ 7 ?M?^ o) = 0, j = 2, 3. 



7,« = 1 



7=1 ?0 = 1 

Coefficients of this system are following (see also formulae (151)): 



(j-l;w;l) _ (j;ti>;2;/9) 
— b al 



093) 



#1 



094) 



(«) 



B^-P™ B[ m -Vi 41) B^-Vi^ 

B W3) _ p(32) s (/34) _ p(42) B <J3j) _ p0'2) 



(j-l;w;m-l) _ (j;w;m;p) 
b aj3 ~ b al 



§(™3) £(m4) 



093) 



-p(32) 
/ a 



094) 



-p(42) 



<$(™j) 

5f > - P^' 2) 



m = 3. 



rp(j-l;wwv) _ rp{j;wv;~il3) 
J a7/3 — 1 al 



B m _ n 



B\ 



093) 



(73) 




- B^ A) 




- B[ lj) 


,(31) 

a 


B (/94) 






_ v m 

1 a 


.(32) 

a 




-p(42) 




T>0'2) 



rp(j-l;wqv) _ j^(lS)rp(j-l;wwv) 



~ Q75/3 



L 7l a<5/3 



p(j'-l;togw;io) 

■ Q7/3 



_ rp(j;ww;/3;i ) 
ayl 



B 0J3) _ p (31) 



(3i ) R (/34) 



(4io) 



5? 4) - Vi A1) 



B^-V\ 



>(/33) 



-p(32) 
/ a 



Jj-l;q;m-l) _ ~(j;w;m;l) 



■ cry 



ay cry 



B\ 

fi(m2) J(m3) fi(mj) 

1 vi' l) vi jl) 

1 V a 32) 



T>(J'2) 



, m = 2,3,j, 





1 


B m 


B 


l;qwu) _ ji(j;wv;jl) _ _ 
a 


1 


-p(31) 


V, 




1 


-,(32) 
/ a 


V, 



rn(j -l ;<?<?") 
' a7<5 



d( 1 <5)'t'(j- 1 ;9' uu ) 



]_ p(3«o) -pO'io) 

i p! 3i) 

1 
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'a/3 



= S 



(j;v;2;aP) 



093) 



B^ B} 



(a3) 



>( Q j) 



'a/3 



0>;m;a/3) 



B|' 



m 



3, J, 



B? 3) - B} 



(73) 



SP 3) - Bi 



(a3) 



il3) 



(«i) 



rp(j~l\VUv) 

1 a-yS/3 



rp(j-l;vuw;io) 
1 a-il3 



rp(j;vw;a/3;io) 
J l7 



093) 



/ 7 



(a3) D(ft') 



p( 1 < 5 )T^(i- 1 ;« OT ) 

^71 J a<5/3 ' 



>(aj) 



,(.j-1;u;to-1) _ g(j>;m;a 



rp(j — l;uvu) rp{ 

J cry — 1 1 



rp(j-l;uuq;i ) _ rp{j\vw;a;i ) 
cry I7 



,j(2m) ^O'm) 



1 B[ 



1 B\ 
1 B{ 



(7i) 



(cy) 



( Q i) 



m = 2, j, 



a7<5 



rt 7l 1 a8 



1 P^ 0) 

1 Bj ai) 



where R^i , B[ 1 ^ and p| tJJ are arbitrary constant parameters. 

Dimensionality of the system (47-51) is defined by the eqs.(45) which involve derivatives with 
respect to 4 independent variables, while other equations involve derivatives with respect to two 
and three variables. Dressing algorithm supplies one arbitrary function of three independent 
variables (which fixes the initial datum for one of the functions p = 1,2), two arbitrary 
functions of two variables (fixing the initial data for v and for one of the functions p = 1,2) 
and single arbitrary function of one variable (fixing the initial datum for u). Namely the system 
(45-51) is written in Introduction, see eqs.(7-13). 

The system (45-51) has the following evident properties: 
d.l) There are infinitely many commuting flows to this system. 

d.2) Equations of this system have differential polynomial structure, also the number of equa- 
tions is rather big. So, one has 24 scalar PDEs in the simplest case Q = 2. 
d.3) Independent variable x\ appears only in the combinations (43), so that it does not increase 
dimensionality of nonlinear PDEs. This happens due to the decoupling of equations for fields 
with different values of the first superscript, which, in turns, is a consequence of the eq.(41). 
However, if one considers another relation ( 40) mixing fields with different values of the first 
superscript, then derivatives with respect to X\ will appear explicitly in nonlinear PDEs which 
will become 5-dimensional. But such system of nonlinear PDEs will have more complicated 
structure. 

d.4) Dimensionality of the nonlinear PDEs is determined by the dimensionality of the linear 
PDE (42) and may be increased without any problem. 

d.5) Reduction to the linearizable system of nonlinear PDEs. The system (45-51) admits re- 
duction v = u = 0, which corresponds to = in the dressing algorithm. In this case the 



>(«) 
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system of PDEs reduces to the single eq.(46) which now reads: 

4 4; « ;4) ^S + E ^VvM + E E ^t^r ii0) = «■ ( 53 ) 

m=l 7=1 io=l 

Eq.(49) is a symmetry of the eq.(53). We may introduce "characteristics variables" d Tct = 
s (4,<j,4)^ _|_ ^^ =i s a' 9 ' m ^dy m , a — 1, ... ,Q which indicate that the dimensionality of the eq.(53) 
is, essentially, min(3,(5). The eq.(53) is " linearizable" because the matrix fields are alge- 
braically expressible in terms of solutions of some linear PDE, see Sec. 3.1. Eq. (53) has been 
written in Introduction, see eq.(14). This equation is also partially integrable in the same sense 
as the system (45-51) does: dressing algorithm provides arbitrary initial condition only for one 
of the matrix fields p = 1,2. 

Proof: The proofs of the items (a)-(c) and (d.l, d.5) are given in Appendix, Sec. 5. The items 
(d.2-d.4) are self-consistent. Relation between solutions of the eq. (53) and solutions of the 
appropriate linear PDE (item d.5) is shown in Sec. 3.1. 

3 Solutions 

In this section all superscripts and subscripts take values from 1 to Q unless different is specified. 
Subscripts do not always denote elements of matrices. We use Greek letters for matrix indexes 
in order to distinguish them from others. 
Solutions of the eqs. (18,19) read: 

$(-»)( A . x ) = j ^ m \\,x)e K(m) ^>^dK, (54) 

C (m) (X;x) = J e Km(qi ' q2 ' x) C^ m \X,q)dq, 

5 

K im \x u x 2 ; x) = (x lXl + k 2 x 2 ) + {P {mj) xi + B^K 2 )xi, 

i=3 

where h = (xi,x 2 ), q = (91,92)- The dressing function must be taken as solution of the 
eq.(15) with k = 2: 

&'\\,H;x) = f (A l/ i;x) + t ) (A lW x) + EW(A^) (55) 

r Q _i 

V { ;\\,^x) = / ^$( sm )(A,x)(x 2 + 9 2 -(x 1 + 9i)A( 2 )) e ^(x^M) x 

CQ m \[i, q)dkdq, 
^\X,/i;x) = j e^^ A(i)xi+Ixi ^(X,fx,>c)dx, 

where \&p^(A, \i\ x) is a particular solution of the eq.(15) and \J/j^(A, fi; x) is a solution of homo- 
geneous equation associated with eq.(15). In addition, we separate matrix function YS s \X,n) 
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independent on x for our convenience, although it might be incorporated in \&j^(A, /i; x). Both 
^^(A, yu; x) and E^^(A, /x) have both diagonal and off-diagonal parts. 
The dressing function G is a solution of the eq.(42): 



G{X;x) = J e^^'^GoiX, x)dx, 

5 

, x 2 ; x) = xixi + k 2 x 2 + 



(56) 



j=3 



It may not be diagonal, otherwise solution space will be poor. 

It is quite standard to assume that the measure dQ(X) has support on an open domain 
V of the A-space, and on a disjoint discrete set of points D = {&i,...,&m}, D D V — 0. 
Correspondingly, we use the following notations for the dressing functions: 



$ (sfc) (A;x) 

C {k) (X;x) 
G(X;x) 



^ sk \X; x) = J <p[ sk) {X, x)e K(k) ^>^dx, XeV, 

<P { : k \x) = ( &{*)e K ' k) ^^dx, X = b n , n = 1, . . . , M, 



(57) 



,(k) 



(A; x) = J e K{k) ^^c (k \q, X)dq, XeV, 



^\ x ) = J e KW to>*rt<$(q)dq, X = b n , n= 1 , . . . , M, 

g(X; x) = J e K(G)( ^'^g {X, x)dx, XeV, 

g n (x) = J e K(G) ^^g nQ {K)dH, X — b n , n = l,...,M, 



U m (X; 



x 



u^ k \X;x), XeV, 
ui^ix), X = b n , n=l,...,M 

and we choose £( S )(A, /j) in the next form: 

f a^(X,fi) = 5 { ; } 5(X-ri, XeV, ^eV 

XeV, fjL — b m , m = l,...,M 
X = b n , n=l,...,M, /j eV 



(58) 



(59) 



(<Tnm)aP = (T nm , X = 6„, fi = b m , 71,771= 1,...,M 



<5 nm , n = l,...,M-l, m = l,...M 
%m, n = M, m = \,...M 

where are scalar parameters. Next, 

^ s) (A,/i;x) = 

4 s) (X, f i;x) = J2 / 0? l) (A,x) J R«(x, 9 ;x)cJ ) (//,g)dx^, X,fieV 
i=i ^ 

i=i ^ 

Vfio^; x ) = ^2 J &(x) R{l \x, T, x)co\v, q)dxdq, X = b n , fi e V 



(60) 



i=l 



r(0 



x)f? w (x, q; x)c£ {q)dxdq, 



X = b n , fi = b r , 
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where n,m = 1, . . . , M, 

flW(x, g; s) = (x 2 + q 2 - (x x + ?1 )^ (2) ) 



(61) 



Remember that the functions $( sm ) are diagonal, i.e. 0( sm ) and 0n"^ are diagonal as well. 

The function y&j^ should be taken in the next form which provides the maximal possible 
richness of the solution space: 



X,fieV 

;L(A;x) = /e^U^ + ^) V ,W 0m(A> x)dx, AeD, ^ = 6 m 

%;no(^ X ) = °> \ = b n , llEV 

^hlm( X ) = > X = b n , fi = b m 



(62) 



where n,m — 1, . . . , M. 

Then eqs. (5) reduce to the following system of equations 

0{7 fc) (A;*) = (yvi 7) (A,i/;x)«^(i/; a ;)di/+ 



(63) 



M 



E^^^ + ^^x))^^) +«S* ) (A;a;), A G T> 

i =1 ' / /3/3 



= I J ^\X,u;x)u {ak) (u;x)du+ 



(64) 



M 



j =1 ' a/3 



(65) 



M \ Q 

Ev&^r^) +e^(#v 

j=l / £0 m=l 



)) 7/3 , n = l,...,M 



= 



^ (zv;x) M ( afc )(z/;x)^+ 



(66) 



M 



Q 



J2^n j (^ ak \x)\ +Y.^{ut k \x))^ U = 1, . . . , M 



3=1 



a/3 



m=l 



for the unknown matrix functions u^ ak \X; x), X G V and u " (x), X = bj, j = 1, . . . , M. 

Following the strategy of [22] , it is enough to provide single linear relation among equations 
(65) and single linear relation among equations (66) in order to satisfy the condition dim ker ^ = 
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I. In turn, this requirement is equivalent to the next two equations: 

M Q M 

E E = 0, £ Af% 3n = 0, (67) 

j = l s=l j = l 

where A^ 71 ^ are some constants. Due to the eq.(67), eqs.(65) and (66) with n = 1, . . . , M — 1 are 
independent equations for the functions u^* k \x), n = 1, . . . , M — 1, where the matrix functions 
u M k \ x ) ma y be taken as arbitrary matrix functions f^ ak \x): u~^ k \x) = f^ ak \x). In order to 
write equations for we involve the condition (41), which reads 

M 

G(A; x) * U^ k \\- x) = g(X; x) * ^(A; x) + ^ ^R^) = (68) 



n=l 



a,k — l,...,Q. 



Thus, the system (65,66), n = 1, . . . , M — 1 and (68) represent the complete system for the 
functions wi^^x), n = 1, . . . , M. 

Once -u( afc ) and u£ k \ n = 1, . . . , M have been found, one constructs the matrix fields -y( mfc ) 
and w nk ' p using the definitions (39): 

/M 
c«(A; x)m^(A; x)dA + ^ cf (x)^ nfe) (x), (69) 

. M 

w^\x) = j g Xp (X;x)u^ k \X;x)dX + J29 jxp (x)u { ; k \x), p=l,2, 



77 



n,i,k = 1,. . .,Q. 

In Sees. 3.1, 3.2, we consider two examples when solutions to the nonlinear PDEs can be con- 
structed explicitly taking into account that the independent variables of the nonlinear equations 
written in Proposition are t, % = 1,2,3, see eqs.(43). 

3.1 Solutions to the system (53) 

Consider the particular case CW = leading to the eq.(53). One has = and ip { ^ nj = 0. 
Thus the system (63-66) reduces to the next one: 

^ fc) (A;x) = (f>a-(A;s)«f^ +u% k \x), XEV (70) 



M 



(3(3 



= E^'^^r^) AGP, (71) 

\J =1 / a/9 

(0i 7fc) (^ = W k \x)U n = l,...,M-l (72) 

= (4 afc) (x)) 7/3 , n=l,...,M-l (73) 

supplemented by the eq.(68) where we take g n = 5 u m without loss of generality. Then eq.(68) 
gives 

u { m\x) = R {ak) - u iak \x), u {ak \x) = g(X; x) * u {ak) (X; x). (74) 
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We need only vr ^ ' with (3 ^ 1 in order to construct the fields , (3 ^ 1. Thus, eqs. (71,73,74) 
yield: 

u% k \X;x) = -{^ m (X;x){r^-u^\x)))^, 0^1. (75) 

The function u {lm \x) can be found as a solution of the following linear algebraic system which 

Q , 

appears after applying / rfA g 77 (A; x)- to (75) and replacing 7 by 7 in the result: 
' =1 Jv 



(76) 



-.(ifc) 



where 



S/T = -(e (0) (A)(i? (lfc) -^ (lfc) )j 7/3 , /3^1, (77) 



Q n 

C = E y rfA ^(A; x) (V$m(A; *)) (78) 



7=1 



V 



(lk;p) 
(lk;p) 



The functions can be found by definition as follows: 



(x) = (g Xp (X;x)*u^ k \X;x)) = (79) 
-(^\x)(R (lk) -u m (x))) , /^1„ p = l,2, 

\ / 7/3 



where 



7/3 

Q 



$ = Ey dA (^(^;^)) . (so) 

7=1 x) 

Now solution q^l of (53) can be found by the definition (44). Formulae of this section suggest 
us to choose M = 1 without loss of generality. 

We see that the fields «/ lfc; p) given by (39) are expressed in terms of p — 0, 1,2. By 
construction, the functions ^ are solutions of the linear PDE due to the fact that g(X; x) is 

solution of the linear PDE (by definition) and functions {^^) Q j{X; x)*j do not depend on the 

variables t, y i: % — 1,2,3 (43) by definition. Thus, the system (53) is C-integrable. 

By construction, functions £^ admit arbitrary dependence on 2 variables yi, i = 1,2. 
Since functions with different values of p introduce the same arbitrary matrix function of 
2 variables, we have Q 2 arbitrary scalar functions of 2 variables and 2Q 2 independent scalar 
fields, which are elements of the matrices q( p \ p = 1,2. Thus, solving, for instance, an initial- 
boundary value problem, only Q 2 scalar fields may be given arbitrary initial conditions. For 
this reason, we consider eq.(53) as a partially integrable system. 

3.2 Degenerate kernel, C w ^ 0. 

The system of linear equations (63-66,68) has a rich manifold of explicit solutions. To construct 
them, we choose, as usual, a degenerate kernel: 

M 

^W) = E^)4?M, (si) 
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where c£j(fi) are diagonal matrix functions. Then 



where 



Q M 

i=i j=i 

Q M 

^pIo(^x) = EE^W( x )4-(^)> 

i=l j=l 



^(A;x) = | ^ ) (A,x) J R«(x,g;x)c;; ) ( ? )rfx ( ig, 



(82) 



In this case, equations (63-66) reduce to the following system of linear equations for the matrix 
functions u^ k \x), u^ k \x): 



MS 1 * 



M Q 



M 



J/3/3 



[ E E rtf (*)*f k) (*) + E W } (*) 

i=l i=l J=l 

(«f fc) (x)) 7 ^, n = l,...,M, 

M Q M 

\ E E ^t lk) (*) + E ^ (*) 



+ 



J=l 1=1 



+ 



3=1 



M Q 



(ti:; ">(x)) ll3 , n=l,...,M, a^P 



M 



j=i i=i j=i 

(uW% p , n=l,...,M-l, 

M Q M 

(^)) 7/3 , n=l,...,M-l, 
k, a, (3, 7 = 1, . . . , Q, n — 1, . . . , M, j = 1, . . . , M, where 



J/3/3 



+ 



a/3 



+ 



(83) 



(84) 



(85) 



(86) 



~(aZfc), 



(87) 



7? 
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and where the given coefficients vf^ 1 , v^ 1 , p£j , p„f , (pn 1 ^ are defined in terms of the dressing 
functions: 



d\ 42(A) (^(A;x) + ^.(A;x)). 



(88) 



v 



Pnj i X ) = ^p;nO;i( X )' 



v 



Eqs.(83) and (84) are obtained applying J dX c ( 2 l(X;x)- to the eqs.(63) and (64) respectively 

v 

from the left. 

Having constructed, from (83-86), the uf k \x), j — 1, . . . , M — 1 and the ^"^(x) in terms 
of u^ k \x), one obtains the functions u( ak \X; x) via the formulae (63) and (64): 



(7*0/ 



*7/9 



Af 



( M Q 

yjyji.<f(A;x )i r>(x) + 

\j=l i=l 



(89) 



09*), 



/3/3 



«^ fc) (A;x) 



( M Q 

EE^ ) (A;-)-r ) W+ 

\j=l »=1 



Af 



E(^'( A - *) + ^(A; s))^*) , A G 2?, a ± (3. 



3=1 



a/3 



(ak) 
M ■ 



Substituting and uf \ j = 1, . . . , M - 1 into (68), one obtains the expressions for 

If g n = S nM , then are defined by the eqs.(74). Equations for the functions appearing 

Q , 

in (74) can be derived applying / dX g^- to the eqs. (89) and replacing 7 by 7 in the 

' =1 Jv 



result: 



m Q 



M 



1=1 i=i j=i 



(90) 



7/3 



_(afc), 



M Q 



M 



) = -EE #° ^t lk) (*) + E (*) . « ^ 



i=l 1=1 



J'=l 



7/3 
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where 

Q 



Z k) ( x ) = E / dX 9m(* x)$ lk \^ x), (91) 
(cfW) = E / dA^n(A;x)(^°(A;x)) , 
(^ X )) 7/3 = E / rfA^ 7l (A;x)(^(A;x) + ^](A;x)) a/3 . 



Now the system (83-86, 74,90) should be considered as the complete linear algebraic system 
for the matrix fields u^ lk \ u^, u^ k \ j = 1, . . . , M, i = 1, . . . , M, k, 1, 7 = 1, . . . , Q where 

uffi = f^ k \ This system is solvable, in general. 

At last, one constructs the fields v ( mk ) ^ w ( nk >p) using (69), and the fields v, u, 

solutions of the nonlinear PDEs (45-51), using the formulae (44). 

Simplest case corresponds to M — M — 1. Then the system (89) yields (using (74) for 

u { %\x) = $ k \\-,x)-(^%\\-x)ur\x)+ (92) 

\i=i 

(^gi(A; x) + ^(A; x))(R^ - ««*>(*))) ^ , A e V 

u [ $\x) = -(j2^l\\;x)u[ mk) (x)+ 
\i=i 

(V'Si(A;x)+^ 1 (A;a;))(it:^)-^(x))) A E V, a^(3. 

/ an 



The eqs. (83,84 ) read (substitute (74) for u { " k) ) 
{^ k \x)) p = [j2^\x)uf k \x) + ^{x){R^ (93) 

i=l 

r Q 1 

= \ J2^ i \x)u { r k \x) + ^\x)(R^-u( ak \x)) +(u { ? lk \x))^, a^f3 

while the eqs. (85,86 ) disappear. At last, the eqs. (90) read: 

*%\*) = iT^)-{fl^\x)uf m \x) + ^ , (94) 

\i=l / 7/ 3 



-(afc)/ 



Q 



' 7/3 



,i=l 



The system of linear equations (93,94) may be solved for {(( aife ) and u^ ak \ Then, the matrix 
functions vS ak ^ will be explicitly given by the eqs. (92). Next, the fields v^ mk \ w^ nk '' p ^ can be 
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constructed using the eqs.(69). Finally, the fields of the nonlinear system can be found by their 
definitions (44). For instance, the explicit solutions in the form of rational functions of expo- 
nents have been constructed. We do not represent them here since they are too cumbersome. 

3.2.1 Richness of the solution space to the system (45-51) 

Let us discuss the richness of solution space, taking into account that the fields of the nonlinear 
equations underlined in Proposition (Sec. 2. 2) are expressed in terms of the scalar fields 

and w£p' P \ i,k, p — 1, ... ,Q, p — 1,2. Using eqs.(63, 64) as the definitions of the function 

u(X; x) and eqs.(69) for the fields w^p'^ and we observe in the small field limit: 

Q M 

~ E^(A;x)*(^^.(A;x) M f ) (x)) i/3 , 0*1, (95) 

7=1 j=l 

Q 

7=1 
7=1 

Q M 

- *(£^(**)«? fc) ^ L 

7=1 j=l 

By construction, the function g has arbitrary dependence on 2 variables y\ and t/2- Functions 
0i and have arbitrary dependence on single variable y\. Finally, iph-,oj do not depend on 

As we have seen in Introduction, the correctly formulated initial-boundary value problem 
involves two arbitrary matrix functions of three variables y^, i = 1,2,3, for the fields 
p = 1,2, three arbitrary matrix functions of two variables yi, % = 1,2, for the fields q( p \ 
p = 1,2, and v and one arbitrary matrix function of single variable y\ for the field u (see 
eqs.(44) for definitions of these fields). However, eqs.(95) show that we have single arbitrary 
matrix function of three variables (formula (95a) after applying J^jj? ■B^' 1 \R~ 1 ) , i 1 g), two 
arbitrary matrix function of two variables (formulae (95b, c)) and one arbitrary matrix function 
of single variables (formula (95d) after applying Ylk p=i '^/? 2 ^(^ -1 )^)- Thus we may not 
supply all necessary initial data in order to completely formulate an initial-boundary value 
problem. For this reason, the system of PDEs underlined in Proposition is considered as a 
partially integrable system. 

4 Conclusions 

We have represented a new type of nonlinear PDEs integrable by a new version of the dressing 
method. This algorithm is based on two principal novelties introduced in [22]: 

1. The nontrivial kernel of the integral operator used in the dressing method: dim ker vj>( sn ) = 
1. This allows one to increase the dimensionality of nonlinear PDEs. In addition, it gen- 
erates arbitrary functions of x in solution space, which allow one to introduce different 



w { ^' p \x) 
v^ k \x) 



(life)/ \ 
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relations among the fields and, consequently, to increase variety of nonlinear PDEs treat- 
able by the dressing method. 

2. The external dressing function G(X; x) which allows one (a) to increase dimensionality 
of solution space and, similarly to the previous item, (b) to increase variety of nonlinear 
PDEs solvable by the dressing method. 

These two novelties seemed out to be extremely promising in development of theory of integrable 
PDEs. This paper is devoted to the investigation of two particular problems outlined in [22]: 

1. Find reductions of nonlinear PDEs derived in [22] which simplify they structure. 

2. Enrich solution space of the constructed nonlinear PDEs (remember that n- dimensional 
nonlinear PDEs derived in [22] have at most (n — 2) dimensional solution space, which is 
not enough for full integrability) . 

We succeeded in simplification of nonlinear PDEs. The systems (45-51) and (53) have differ- 
ential polynomial form, although the number of equations is rather big. We also have enriched 
solution space. Thus, solution space of the constructed 4-dimensional PDEs is parametrized 
by one arbitrary matrix function of 3 independent variables, by two arbitrary matrix function 
of 2 independent variables and by one arbitrary matrix function of single independent variable. 
However, we need one more arbitrary matrix function of 3 independent variables and one more 
arbitrary matrix function of two independent variables in order to achieve full integrability. In 
other words, we achieve only partial integrability of the derived nonlinear PDEs. 

A natural generalization of this algorithm is an increase of dimensionality of the kernel of 
the integral operator: 

dimker ¥ sn) > 1. (96) 

The study of associated (partially) integrable equations is postponed to future investigation. 

Author thanks Prof. P.M.Santini for useful discussions and referee for important comments. 
The work was supported by INTAS Young Scientists Fellowship Nr. 04-83-2983, RFBR grants 
04-01-00508, 06-01-90840, 06-01-92053 and grant Ns. 7550.2006.2. 

5 Appendix: proof of Proposition, items (a-c) and (d.l, 
d.5) 

As we have seen in Sec. 2. 2, the eq.(38) leads to the linear relations between any two solutions 
jj{nk)^ ^(ji-M)^ £(j3;nfc)^ j > 2 and L^ 2 ' nk \ j > 3 of the homogeneous integral equation (31). 
In orther words, one has the following system of independent linear equations for the spectral 
functions U^: 

L { S nk) = £ L^F%™ =► (97) 

il,7=l 

f) TT( nk )j_ft rr( nfc )K?( fc J» f) jMk) rj(kj) , 
11,7=1 H,7=l 
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L ( S' lk) = E L%MF%™ =► (98) 

il,7=l 

Q 

°xiU a p A n + 2-*t a 7 7/3 ^7 / D /3 _ 

il,7=l 

E ugw™$?* k \ j > 2. 

il,7=l 

^s ;/fc) = e 4 2 7 i;Hi) ^r ife) =► (") 

il,7=l 

Q 

Tj(nk)„{j;n) _ ST~^ n(nii) »(2;ra) ft{j;hk) 
U al3 °/3 — 2.^1 a ~f 7 ^7/3 
il,7=l 

where F^ lk , F^ lk and F^ lfc are scalar functions of a; to be defined. This set of equations 
can be taken as a basis. It is simple to check that any other linear equation for the spectral 
functions f/( nfc ) is a linear combination of (97-99). Remark, that if n — /?, then LHS of the eq. 
(98) disappears yielding algebraic relations among elements of matrix spectral functions: 

E U^B^P\f k) = 0, (100) 

il,7=l 

which must be involved into consideration as compatible constraints. We consider equations 
(97,98) and (100) as an overdetermined system of linear equations for the spectral functions 
jj(nk) di sre g ar di n g equation (99). It will be shown in the Sec. (5.1) that eq.(99) follows from the 
eq.(100). 

Now we are going (a) to derive nonlinear equations for the fields (39) and (b) to express 
functions F(x), F(x) and F(x) in terms of the fields (39). 

Nonlinear equations for the fields v ( mk ) )V (. mk > 1 ) result from the eqs. (97,98,100) after applying 

Q 

E/ ^&a* t° them and using (19) for Cx}, k > 2 (we replace a by a in the result): 



a=l 



Q 

(nik;l) p(v;jl;nik) , Wik;2) p{v;j2;ik) \ ^ „,(mii)„,( ni i fc ) ( r>(hj) u( k i)\ — 
a/3 °a/3 V a/3 °a/3 U a~/ u 7 /3 y n -y n /3 J ~ 



V 

il,7=l 

Q 

E v^B^F^ k \ j>3, 

il,7=l 
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d,/:; k) - d xl v% k) A^ + (102) 

Q 

} (nik;l) ^( v -jl-ni) + v (nik-,2) g( v -j2-,i) + ^ ^(nji^Mlgfii j)^ gO 2 ^) _ 

il,7=l 



il,7=l 



£ v^BMpff^ = 0, (103) 

il,7=l 



where 



5 (j, y -l;nifc) = _ ^;n) = ^) _ ^) + ^fe) ^(2) _ 3(^(2) ^ (1Q4) 

e(v;j2;ik) _ R (fcj) _ R (jj) 
o(i>;21;ni) _ A (2) o(>;22;i) _ _ 1 

Sgtfi;™) = A W + gfeia) = A W _ A 0') + Bfe)^ 2 ), = -B^\ j > 2. 

In order to fix functions !F^' Mk \x) we introduce the external dressing matrix function 
G(\; x) and the associated fields 



it; 



= G* U (nk \ w {nk]p) = G Xp * U (nk \ w {nk]ps) = G XpXs * U (nk \ . . . , w {nk]ps) = w [nk ' sp) (105) 



(see eq.(39) in Proposition ). Applying Yla=i G&a* to the eqs. (97,98, 100) one gets the first set 
of equations for these fields (we replace a by a in the result): 



Q 

_„.,(nk;l)r>(.kj;n) (nk;2) R (kj) (nil) (niifc) / R (i u ) E>Oi)\ _. 

"V °/3 ~+~ ^ -°/3 ^ " ; a7 U 7 /3 ^-° 7 J ~ 



ii,7=l 



w™B^F$ lk) . 



ii,7=l 



(S^™ - 3^ ^) - (107) 



J2 w^B^F^ k \ j>2, 

il,7=l 



EEW^f^O. (108) 

«i=l 7=1 
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Applying Yla=i(Gaa)x p * to the eqs. (97,98, 100) and replacing a by a in the result, one gets 
the second set of equations for the fields (105): 

d Xj w% h ' p) + d Xl w% k -' p) Bf j -' n) - d X2 w% k -' p) Bf j) - w% k ^ (109) 

Q 

(nk;pl) K (kj;n) <nk; P 2) R {kj) \ ^ (ni i; p) (niifc) / p (i x j) f>(kj)\ 

W a/3 D /3 ~+~ W a/3 n /3 ~+~ W a-y U -y/3 ^-°7 " 



H,7=l 



£ ^CniiSP)^)^*). 
il,7=l 

(^^ M -^^^)- (110) 

Q 

..,{nk;pj) , (nk;pl) A (j) 7 „(mi;p) fan*) n(iij) \ »(2;n) _ 

ii,7=l 

Q 

^ ^isp^F^, j>2, 

U,7=l 

£ w^E^F^ k) = 0. (Ill) 

il,7=l 

We refer to the eqs. (106-111) as equations for the fields u>( nfc;p ) and w^ nk ' ps \ also the nonlinear 
parts of these equations involve fields i)( mfc ) as well. 

Q 

The set of nonlinear equations may be continued applying ^(Gfi tt ) Vs * to the eqs. (97,98, 100), 

a=l 

but the system of nonlinear PDEs constructed in this way may not be completed since (a) its 
solution space has arbitrary functions of all variables x, (see functions in the eq.(38)) and 
(b) the dressing function G is arbitrary, i.e. all derivatives of G with respect to x p , p — 1, 2, . . . 
are independent. This ends the proof of the item (a) of Proposition, i.e. we have derived a set 
of nonlinear PDEs (101-103), (106-111) for the fields (39) which does not represent a complete 
system of PDEs. 

Now we pay attention to the completeness of the derived nonlinear system. 
In order to fix arbitrary functions in the solution space we must introduce an extra largely 
arbitrary relation among the fields 

T rs (all fields ) = 0, r,s = l,...,Q, (112) 

where T rs are the Q x Q matrices. (See the item (b) of Proposition). For instance, in order 
to obtain a system of nonlinear PDEs having differential polynomial form, we choose the eq. 
(112) in the form (41) (see the item (c) in Proposition ). 

We use the eqs. (106,107) to express F^ nk \ Pfy 11 ^ in terms of the fields (39). In view of 
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the eq.(41), the eqs.(106,107) read: 

Q 

-™% k ' d) ~ w^Bfri + w™*>B™ + R [ S l)v ^ k) { B i lj) - B t j) ) = ( 113 ) 

il,7=l 



E R^B? n) F$ hk \ J > 3, 

ii,7=l 

( - ^ + ^ n ; ;i) 4 j) + e R^t )B ? j) ) B f n) = 

il,7=l 

E B^B^F^ k \ j>2, n^p. 

il,7=l 

The matrices must provide unique solvability of (113) with respect to F^ 11 ^ and F^ nk \ 
which become linear functions of the fields v and w: 

Q 

ip(j;lk) _ q(j;Uli2k) (ill2fe) i hliN 

U,l2,7=l 

E c* -1 )^ ( - w< if ,j) - ^f' 1)B f' ,il] + ^f ,2)B t j) ) > ^ 3 ' 

il,7=l 

Q 

a/3 — 2.^1 "7/3 7/3 " + " 
U,i2,7=l 

E ( - ™§ M + ^ fc;1) ^>? n) , J > 2, 



il,7=l 



where 



= E(^ 1 )S 1 1 ) ^717 2) (^ 2J) " *f\ ( 115 ) 



71=1 



71=1 

and the operator R" 1 is the inverse of the operator R^B^ 2 ' n \ i.e. 

Q 

E {R- 1 )^^^ = <^ } <W or (116) 

n,7i=l 



E < 7 i l) ^(^" 1 )7S = ^ n0 ^a- 
U,7l =1 

The eq.(lOO) may be transformed as follows. Substituting F^ lfe ) into (100) one gets 

E U^B^\R-^B^£^ = (117) 

U,i2,7l>72=l 
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where 

Q 

c(i2j) _ \^ D(i 2 i 3 ) ofei), ,0^3*0 _ (i 2 fey) , (i 2 fe;l) /-, , ox 

fc 72fe - ^7273 ^73 V 73/3 ^72/3 + ^72/3 ^2' 

23,73 = 1 

i 2 = l,...,Q, j = 2, ...,Q + 1, fc,7 2 = 

Assuming large arbitrariness of v^ 3k ^ and, consequently, of we conclude that coefficients 

ahead of £^ 2 ^ are zeros: 

E U^B^\RT^B^ = 0, (119) 

»1.71 = 1 

We take i 2 ^ (3 in (119) since LHS of this equation is identical to zero if i 2 — (3 due to the 
equality B^'® = 0. Thus we have Q 3 (Q — 1) equations and the same number of the elements of 

the spectral functions U^^\ 71 7^ (3 in (119). Remark that eqs. (119) are linearly dependent. 
Namely, they admit Q linear algebraic relations for each pair (a,/3). To show this, one applies 

E^iE^i'^S?' & = to the eq.(119). This yields identity due to the eq.(116). 

One concludes that the system (119) with fixed a and (3 consists of Q(Q — 1) — Q equations 
and consequently Q(Q — 1) — [Q(Q — 1) — Q\ = Q elements of the spectral functions Jjf^\ 
71 7^ (3, may be independent for each pair (a, (3). If (3 ^ Q, then ujfg 1 ^, i\ = 1, . . . , Q may be 

taken as independent functions in the set U^i \ H, 7i = 1, • • • , Q, 7i 7^ P- 

Q Q 

Applying C^l* and ^^(Ga a ) Xp * to the eq.(119) and replacing a by a in the result one 
gets the following algebraic relations among the fields v^™ 1 ^ 1 and Wajl' P ^' 

E v^B^\tC^B^=0, (120) 

*1.71 = 1 
71^/3 

Q 

E ^^^fKR-^fBf^ = 0, (121) 

n,7i=i 

i,i 2 ,a,(3,l2 = 1,...,Q, i 2 + (3. 

The last step which must be done to complete the system of nonlinear PDEs is introducing 
additional PDEs for the dressing functions G establishing relations among derivatives G Xp and, 
consequently, among the fields w^' ,p \ The simplest case leading to 4-dimensional nonlinear 
PDEs is following: 

d Xk G = V {kl) d Xl G + V {k2) d X2 G, k > 2. (122) 

Then 

w (ij;k) = V kl w (ij;l) + V k2 w (ij;2)^ k>2 ^ ^3) 
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Not all equations of the system (101-103), (106-111) are independent. Namely, although 
the linear equations (97) can be written for any spectral function U^ nj \ it is not necessary to 
use all these equations because the eqs.(97) with n ^ (3 are combinations of the eqs.(98) by 
construction. Thus, we need only those equations (97) which correspond to n — (3, while the 
equations for U^J\ n ^ (3, are represented by the system (98), which is simpler. Since the 
eq.(97) generates nonlinear equations (101), (106) and (109), one needs only those of them 
which correspond to n — (3, while the rest equations of the nonlinear system will be generated 
by the eq.(98), see eqs.(102,107, 110) 

Now we write down the complete system of nonlinear PDEs for the fields (39) using 
eqs.(120,121) and (116) (see eqs.(124-137)). Thus, eq.(101), n = /3, yields: 

ti a/3 ■- OxjVap + Xl V a/3 L5 p - X2 V a/3 Hp + 

mk;l) q (v;jl;/3ik) (Pik;2) Mv;j2-ik) 
U af3 °a/3 V a/3 °a/3 



where 



E~,(Pih)( Q (v;j;hi3k) \^ \^ (l3k;i ) q (w-,i -,j;hk)\ _ n . . „ 

V «li \ V 0P D 7l/3 + W 72/3 ^7172/3 / ' 3 ~ ' 

H.71=l *3 = 1 72 = 1*0 = 1 

s tehi li3 k) = ^(m 3 /3) _ 5 (i3h) 5ii ^ S$*' iak \ (125) 
Q 

gfrisn) = J2(R-%^R^\ J > 3, (126) 

72=1 

S^f lk) = B^\R-%^S^ r ' k \ (127) 

q (w;2;j;hk) _ »(2;/3) / 6-l\ (ii/3) q{w;2;j;k) / 19 o^ 
D 7l72/3 ~ °7l y U ^7i72°72/3 ' 

5 (j»;l;i;fc) = _^1) _ ^) = _pO'D _ B f j) Af + A^, (129) 

S yM = -vW+Bf\ j>3. 

Writing the expression (125) we used eq.(120) and eq.(130b) from the following list of identities: 

Q Q 

U^B^ n \R~ 1 )^R^ 7l ^n, => (130) 

U,7i72=l U,7=l 
Q Q 

U,7i72=l U,7=l 
Q Q 

H,7,72=l «l,7=l 

n,z,a,7i = 1,...,Q, p = l,2, ii^n. 

Q 

To prove equality (130a) one must apply ^ ■B^ a) (Rr 1 )^ to (130) and use (119) to simplify 

71 :*3 = 1 

the RHS of the resulting equation. Eqs.(130b,c) follows from the eq. (130a) after applying 
Q Q 

Cj?1* and y~^(G f a a ) Xp *, p = 1, 2 to the eq.(130a) and replacing a by a in the result. 

a=l a=l 
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fp(j;nik) 



The eq. (102) yields 

{d*^ k) -d Xl v$ k) AP+ (131) 

(mfc;l) q(v;jl;ni) , (nik;2) q(v;j2;i)\ »(2;ra) 

Q Q Q 2 

E Jnih) ( \ (ni 2 fe) q(vy;mii 2 ) \ \^ (nk;i ) q{w;i ;j;nh)\ ■ > 9 / o 



ii,7i=l 12=1 72=1 i =l 



r,(w;l;j;rai) _ »(2;n) / p-l\ (iin) »(2;n) c(w;l;j;n) ( i qq\ 

^7172/3 ~ °7i ^ ^7172 D /3 D 72 ' 

^S' ;mi) = B^ n \R~%^Bf' n) S^ n \ (134) 



where 



^;l;2;n) = A (2) ? ^;2;2;n) = _^ ^35) 
cO;l;j» _ _-T?0'l) I A(j) C(w,2;j;n) _ _<r>(j2) „' > o 

Writing the expression (132) we used the eqs.(120) and (130b). 
The eqs.(109) with n — (3 gives 



2 

E„,,(/3fc;P*o) c{w;i o ;j;k) K (kj]0) 
W a/3 °al3 °/3 ~+~ 



i =l 

Q Q 



E,„(/3»i;p) f „ (0*2*0 q( u >r,hi2k) , \ (/3fc;i ) c(w;ioy';iifc)\ _ n 

W «7i ^ 2^ U /3/3 ^71/3 + ^72/3 ^7i72/3 J ~ U ' 



U,7l=l *2=1 72=1 «o=l 

Finally, the eq.(llO) reads: 

;= ^(J*) _ d^^) + £ u,^)^"*^)) 4 2;n) + (137) 

io=l 

E ( E <T fe) s 7 if ni2) + E E ^3 ;i0) ^S j;mi) ) , j > 2, « ^ (3. 

U,7l=l *2=1 72=1*0=1 

Deriving the coefficients in the nonlinear parts of the eqs.(136, 137) we used the eqs. (121) and 
(130c). 

Now we consider such combinations of the eqs. (124), (131) which do not involve the fields 
v (mk;j)^ j _ i 2 ; and such combinations of the eqs. (136) and (137) which do not involve the 
fields w^ nk ' ,ino \ p — i 2. There are several types of such combinations depending on the 
relations among the coefficients ahead of the fields v^ mk '^ in the eqs. (124), (131) and ahead of 
the fields w ( nk ' pi ^ in the eqs.(136), (137). Remember that the parameters S { ^ jl '' ili2h \ S { ^ j2 '' ili2 \ 
Sat ni2 \ S { :f' h \ S^' j *\ S^ k) and B ( : r ' h) are given by the eqs.(104,129,135,27). 
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1. Let n = = a, i = k. Then sj£ liakk) = Sj£ 2ikk) = 0. The eq. (124) reads: 

rp(jl;akk) — ft „,(akk) , ft (akk) K>(kj;a) _ ft Aakk) p>(kj) , 
^aa - — u Xj u aa ' u xi u aa u a u x* u nn -D« ^ 



^2 "aa -^a 

Q 2 



(138) 



Q Q Q 2 

E (afeii) / (ai 2 fc) o(v;j;hi2k) , (afc;i ) oO;«o;j;nfc) ) _ n 

a 7i V / / «a ^7ia ~ / , , / . , "^72 a 7172" y 



il,7l=l 



«2=1 



72 = 1 i =l 



2. Let /3 = a in the eq.(124) and n = a in the eq. (131). Then S& y ' 1;aifc) = -Ai 2) S^' 2;ifc) , 
S^ jl ' a ^ = —A^S^p 2 ' 1 ^- So that the following combinations of the eqs. (124) and (131) have 
no fields v%% k '^: 



£j(j2;aik) ._ 



rp(3;aik) rp(j;aik) 



(v;32;ik) ~(v;j2;ik) 
*-^aa 



E 



(j3;aik) 
a/3 



s, 



^{v;22;i) Mv;j2;i) 
Oct Oa 



, j>4, 
, J>3, 



(139) 



which have the following explicit forms: 



j^(j2;aik) 



8 vwM dxjV g*) + £ a g>^*) 5iBmV («*) _ (140) 

m=l 

Q Q 

E (aih) (ahk)rp(j;vv;ihk) , \ A (aiii) (012k) c{iii20c)rp{j;vv;ii2k) , 
«a "aa a / ,, 07 aa ^7 a ~r 



i l> i 2'T = 1 



where 



il, 71. 72=1 10=1 
71 



F (j3;aik) ~(j;v;j;i)a „M0 , \^ 5 0>;m;i)o (< 



aik) 
■fi 



m=l 



Q 



Q 



^ v (aiii) v ( a ii k )f<(3-,vv-,ih) _|_ ^ v (aih) v (.ad 2 k) gfa^a) rp(j;vv;ii 2 ) _|_ 



i i= i 



H,»2.7 =1 
«2? ii '77 £ a 



Q 2 

E E ^^S^^H^" 1 )^^*^ = °> 



i 1.71.72 = 1 10=1 

J > 3, a^P, 



(j;v;l;ik) 



(fc3;a) 



(kj;a) 



R (fc3) _ R (i3) R (fcj) _ R («) 
-Da JJa -Da -Da 



(j;u;2;ifc) 



(j;v;m;ik) 



B. 



(fc3) 



5, 



R (fe3) _ R (i3) R (fej) 
-Da -Da Da 



5, 



j(3m) 



R (fc3) _ R (i3) R (fej) 
-Da -Da -Da 



5, 



, m = 3,j, 



(141) 



(142) 
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rp(j;vv;iiik) _ 



R (fc3) 
J->a 

R (fc3) 



d(ii3) 

XJq; 
R («) 



r>( fc i) o(uj) 

-Da -Da 



rp(j\vv;iik) rji[j-vv;kiik) rp[j^vv;ikk) r\ 



rri(j;vw;ik;l) 
07 



rp(j;vw;ik;2) 
07 



£(fc3;a) + p^31) 



R (fc3) 



R (i3) 



R (*j) 
-Da 



d(«) 
-Da 



s (fcS) _ p (32) 

R (fc3) R (t3) 
-Da -Da 



-Da -Da 



rp(j;vw;ii;l) _ n 
^07 ~" u - 



rp(j;vw;ii;2) 
07 



0, 



A a 

ji(j;vw;i;l) 
07 



4(2) jC?) 
^T-a ^la 

1 fl^ 
1 

1 



j(j>;m;i) 



£(2m) £(jm) 

1 



4 (2) jO") _ T)(J1) 
1 flj«> 



ji(j;vv;ii) _ g 



- 07 



, m = 2,j, 



1 > 
1 



(143) 



3. As for other values of n, i, k, a, (3, the following combinations of eqs.(124) and (131) have 
no fields v^ k '^: 



E 



a/3 



E 



(3;(3ik) 
a/3 



E 



(4;/9tfc) 



g(v;31;/3ik) g(v;41;f3ik) g 



Q(v;32;ik) 



a/3 
;;41: 

v;A2;ik) 



E 



(j;f3ik) 



af3 
(v;jl;/3ik) 



v;j2;ik) 



'aj3 



, i>5, 



E 



(jb;nik) _ 



a/3 



E 



(2;nifc) 
a/3 



p(3;nife) j-,{j;nik) 
^ap ^a/3 



r>(u;21;ni) o(«;31;ni) r>(u;j'l;m) 
Oq, O a >5q 



Q(f;22;i) 
Oa 



qO;32;j) 
Oa 



c(«;i2;i) 
Oa 



, J>4, 
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These expressions have the following explicit forms: 

4 



E 



(j4;(3ik) ._ (j;v;j;ik) o (/3ifc) (j;v;m;ik) r, (/3ifc) 

a/3 •— 6 a/3 U Xj U a(3 Z-^i a P °^ U a/3 

m=l 



where 



nXPih) A/3iik)rp(j;vv;iiik) 
V af3 U Pf3 1 a/3 



U a~/ U f3l3 C j 1 al3 



E 

»1.T1 = 1 



*1 '*2 j7— 1 



<1,71.72 =1 «0 = 1 



^a/3 •— S a % V a/3 + a x "> 1 



(nik) 
a(3 



m=l 



Q 



Q 



E (niii) (mik)rp(j;vv;iiin) , \ r ,,( ni 2fc) c(iii2n)np(j-,vv-,ii2n) , 

V y n/3 J a a 7 "Z 3 S i « ~ l ~ 



il=l 



»1><2'7=1 
7#" 



Q 



E E^^S"^?"^^ 1 )^^ 



vw\in;i 
072 



0, j>4, 



<1,71.72 = 1 10 = 1 
71 



'a/3 



„0>;2;ifc) 

'a/3 



»(fe3;/3) e(fe4;/3) 
°/3 °/3 
o(«;31;/3ifc) ~(i;;41;/3ifc) q(v;jl;/3ik) 
^a/3 ^a/3 <-a/3 

Q (v;32;ik) Q {v;A2;ik) Q (v;j2;ik) 
^a/3 ^a/3 ^a/3 

R (fc3) R (fc4) R (fcj) 

-°/3 n /3 n /3 

q(v;31;/3ik) ~(v;4l;f3ik) q(v;jl;/3ik) 
^aB ^a/3 ^a(3 

q(v;32;ik) ~(v;42;ik) q(v;j2;ik) 
D a/3 °af3 °af3 



Jj;v;m;ik) 
5 a/3 



£(lm) 

c O;31;/3ifc) 
^a/3 
<-.(«;32;jfc) 
^a/3 



£(2m) 

c (i);41;/3ifc) 
^a/3 
c (-u;42;jA:) 
^a/3 



£(3m) 

(v;jl;/3ik) 
a/3 
c(v,j2;ik) 
^a/3 



5. 



2, 3, j, s 



(j;v;m;ik) 
a/3 



>f3a 



p(j\vv;iiik) 
[ a/3 



<-.(u;32;iifc) 

Q (v;31;f3ik) 
^a/3 
Q(v;32;ik) 
^a/3 



<j(u;42;iifc) q(v;j2;iik) 

^00 ^(3/3 

<-.(« ;41;/3ifc) q{v;jl;(3ik) 

a/3 a/3 



,(?;;42;iA;) 
^a/3 



Jv;j2;ik) 
5 a/3 



rp(j;vv;ikk) _ n 
J a/3 — U > 



■.(«j;io;3;fc) <-.(uj;io;4;fc) q(w;i ;j;k) 



rp{j;vw;ik;i ) 
J a7/3 



^7/3 
c (i);31;/3jfc) 

<^a/3 
<->(i;;32;jfe) 

^a/3 



^7/3 
c (i;;41;/3ife) 

^a/3 
c (i>;42;ifc) 

^a/3 



^7/3 
c (t);jl;/3jfc) 

^a/3 
c (v;j2;j/c) 

D a/3 
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Ah) 

Sin 


Ah) 

Sin 


Ah) 

Sin 


~(j;v;l;ni) _ 




£(v;31;ni) 
O a 


Q(v;jl;ni) 
O a 




g(v;22;i) 
O a 


£(v;32;i) 
O a 


<;(v;j2;i) 
O a 



S( lm ) §( 2m ) 



(J(3m) 



~(j;v;m;ni) 



Q(t);21;ni) X(v;31;ni) &(v;jl;ni) 

Oa Oa Oa 

q(v;22;i) £(v;32;i) A(v;j2;i) 

Oa Oa Oa 



m 



3, 4, j 



~(«;22;«i) t ,(i);32;i 1 ) ~(«;jr'2;ji) 
O n O n O n 



On 



(v;21;ni) g(v;31;ni) Q(y;jl;ni) 



So 



£(,v;22;i) ^(v;32;i) ^{v;j2;i) 
Oa Oa Oa 



rp(j;vw;in;i ) _ 
«7 



q(w;i ;2;n) q(w;i ;3;n) q(w;i ;j;n) 

O'y O'y O'y 



(v;21;ni) 
:(v;22;i) 



;(u;31;ni) 
O a 

(?(v,32;i) 
O a 



*(v;jl;ni) 
'a 

:(v,j2;i) 



As we shall see in the end of this section, the eqs.(144,145) do not appear in the final complete 
system of nonlinear PDEs. For this reason, we leave coefficients T, T, s and s in general form. 
4. As for the eqs. (136) and (137), the following combinations should be considered: 



E 



(j6;(3kp) _ 
a/3 



E, 



(j7;nkp) _ 
a/3 



p (3;/3fc;p) 

c (w);l;3;fc) 

<-.(mi;2;3;A;) 

°a/3 

fi(2;nfc;p) 

a(3 
<r(w;l;2;n) 
O a 

~0;2;2;n) 
O a 



p (4;/3fc;p) 

c O;l;4;fc) 

Mw;2;4;k) 
^af3 
p(3;nfc;p) 

£(w;l;3;n) 
O a 

~0;2;3;n) 
O a 



^a/3 

c (io;l;j;fc) 
D a/3 

Q {w;2;j;k) 
°af3 
£p(j;nk;p) 

<r(w;l;j;n) 
O a 

£(w;2;j;n) 
O a 



i>5, p=l,2 



j>4, p=l,2, 
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These expressions have the following explicit forms: 



where 



rp(j%;f3kp) 
^a(3 



E 



(j7;nkp) 

a/3 



Q 



m=l 



i (f3h;p) nt (f3hk) r p(j;wv;iik) 
-a/3 



/3 u /3/3 J 



i 1= i 



il>i 2 >7 =1 



,(nfc;p) 
*/3 



e e ^^^^rn^ 1 )^^*^ = o, j > s, 

n,7i,72 =1 io=l 

71^/3,725^ 

3 

rra=l 

Q 

E w (nu;p) u ( ni i fc )j'0>f;nn) _|_ 

ii=i 
Q 



H. < 2.7= 1 
75^» 



E E <rM:? 0) #? n) (^ 1 )S5 ,8in;i9) = o, 

il,71.7 2 = l lo=l 
71^".72^« 

J>4, 



'a/3 



Sr,w;2;k) 
'a/3 



13 



LJ 

(km 
& 



(fc3;/3) 



p (31) i3 ( fc 4;/3) +pr ; g 



»(fc4;/3) 
°0 



(41) „(fci;/3) 



R (fc3) _ —(32) 
R (fc3) 



S? 4) - Pi 42) 



/3 



B 



(kj) 



B 



(k4) 



5 



R (fc3) _ —(32) 
/3 ' a 



R (fc4) _ —(42) 
/3 ' a 



5 



-pO'2) 



^(j;w;m;k) 
J af3 



£(m3) 



£(m4) 



ft(mj) 



g(k3;f3) + p(31) g(fc4;/J) 



(41) 



s (fc3) _ ^(32) 



^fc4) _ ^(42) 



K (kj\P) 



B (kj) _ pt?2) 



3,4,j 



T 



(j;wv;hk) 



a/3 



R (fc3) R (ii3) 



B (3 - B 3 



B (km + V m B 



B 



(k3) 



-p(32) 



5 



(fc4) 



T)(42) 
' a 



(km 



B 



(kj) 
P 



- 

I a 



T. 



(j;wv;kk) 



a/3 
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p(j;ww;k;l) 
-Q7/3 



p(j;ww;k;2) 



~(j;w;l;n) 



(fc3;/3) 
3 

(fc3;/3) 



p(31) 



7 

(31) 



B (fc3) _ p(32) 



+ p (41) 



K>(kj;t3) 
°3 



Cji) 



5 



(fc3) 



-p(32) 
/ 7 



5 



(fc4) 

/3 



-p(42) 
/ 7 



B 



(kj) 



I 7 



B m P) + p (4D 



R (fc3) _ ^ 



(32) 



R (fe4) 



(42) 



,(2) 
4(2) 



~(j;w;m;n) _ _ 



rp(j ;wv;iin) 



£(m2) 



J (3) 

4 (3) -p(31) 

T>(32) 
/ a 

£(m3) 



J (J) 

-pO' 1 ) 

T>0'2) 
/ a 



.4 



(2) 



4 (3) _ -p^ 



(31) 



£(™i) 



(32) 



1 

4(2) 

1 



r(*i3) 
Dn 

4(3) t>(31) 
■r^n r a 

<p(32) 
/ a 



7>, 

r(*U') 



(J2) 



2, 3, j, 



4O) _ -pu 

p0'2) 
/ a 



Oi) 



^0';«"";fe;i) _ n 

J aa/3 — U ' 



rp(j;ww;k;2) _ n 
J aa/3 — U 



rp(j;ww;n;l) 
a-y 



rp(j;ww;n;2) 
07 



4(2) 
,(2) 



1 



4^ 



(3) 



p(31) 



-p(32) 



<p(32) 
/ 7 



(2) ,(3) 



4 l ; 4 



<p(31) 



p(32) 
/ a 



4(i) 

^(i) _ p(il) 



p(il) 
/ 7 



X>(J'2) 



i0) 



(Jl) 



T>0'2) 



?fi(j;tuiu;n;l) _ n 



nn{j;ww;n;2) _ n. 



It is remarkable that the system of nonlinear PDEs (138,140,141,144, 145,146,147) has the 
complete subsystem of nonlinear PDEs. Let us show this. First of all, we remark, that the 
fields v^aa^ appear only in the eq.(138). Moreover, the fields v^ k \ n / a do not appear 
in the eqs.(140, 141, 146, 147). Thus namely these equations compose the complete system 



of nonlinear PDEs for the fields v^ k ^ and This system is compatible with the system 



(nk) 



J aj3 



J a/3 



of algebraic relations (120,121). Eqs.(144, 145) must be considered as linear PDEs for the 
functions v£* k \ n ^ a, with coefficients depending on the fields v^g® and w^t\ In addition, 



J a/3 
(aik) 



the fields and w^ g k ' p ^ with different values of the first superscript appear in the decoupled 



J a/3 



subsystems. For this reason, we fix the value of the first superscript in the fields to be 1. For 

may be taken as independent fields, while 



(life) (life) (lfc;p) 



IQ 
(lk;p) 



and w { ^ p) 



(lk;p) 



this choice, the fields 

the other fields and w^' p \ (3 7^ 1,Q, are expressed through v^q^ and w ( ^q P> owing to 

the eqs. (120,121). However we will use the dependent variables u a g and q^l defined by the 



eqs.(44) instead of Vy£^ and w^q P> respectively in order to simplify the final nonlinear PDEs. 



fife;?) 



'1Q 



J aQ 



This choice allows us to remove algebraic relations (120,121) from the final complete system of 
nonlinear PDEs. Finally, we observe that, for our choice of the first superscript in the fields, 



parameters A± \ j > 2, appear in the combinations Af — Va L> ■ For this reason we may put 



(3) 



>Cji) 
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= 0, j > 2 without loss of generality. To simplify formulae, we assume that the matrix 
A^ satisfies the following two requirements: 

1. 4 2) = -1, (149) 

2. Af^Af, i^j, i,j 1 Q. 

Q 

All in all, applying -B^ 1 ^ -1 )^ to the eq.(141) with a = 1 and to the eq.(147) with 

/3,fc=l 

n — 1, putting a = 1 into the eq.(140) and [3=1 into the eq. (146) we end up with the system 
of nonlinear PDEs (45-51) for the fields (44) (we write the symbolic set of equations in the 
same order as they are written in the system (45-51)): 

£(56;w £( B7ilA0 = J2 E^B^\R-%f, (150) 

7,fc=l 
7,fc=l 

^(47; W = ^ £ C47;l* 1 ,) B (2;l) ( ^-l ) (*l) j 
7,fc=l 

/3,fc=l 

where we have redefined constant coefficients in order to improve the structure of PDEs: 

(j-l;w;m-l) _ (j;w;m;/3) _ 9 o a ■ rp{j -l;wwv) _ rpfrwvnP) fici^ 

S af3 — S al > 171 — Z ?<J) 4 !J) i a 7 /3 ~ J al ) l i0i J 

^(j-ljtugwjio) _ T (j;ww;/3;i ) • \ r 
■ 1 a7 / 9 — J a7l ' J — ° 

Jj-l;q;m-l) _ ~(j;w;m;l) _ 9 q • T->(j-l;^u) _ ^(j'siuv^l) 

7 i O'-i;9??;io) _ 7 : '(i;"'w;i;j ) ,,• > 4 

- 1 «7 - 1 07 ) J — ^> 

(j"-l;u;m-l) _ (j ;v;m;a/3) _ n „• rp(j-l;vvv) _ rp{j;vv;a-yl3) 

■ b a/3 ~ *1 > "I' — °i Ji J a7/ 3 — J l ) 

rpCj-ljtmwjio) _ rp(j;vw;al3;i ) ■ a 
1 a. 1 fi — J l7 ' ^ — ' 

0-l;u;m-l) _ ~0>;m;a) _ 9 • T (j-l;ut)u) _ f>(j;™;a 7 ) 

*a — *1 ) " L — Z )J) J a 7 — J l i 

T->(.j'-l;uug;io) _ rfi(j;w«;a;io) „• -> o 
J «7 — 17 ' J — °" 

This ends the proof of the item (c) of Proposition. The eq. (53) of the item (d.5) follows from 
the eq. (46) if one puts v n k = u n k = 0. Existence of the infinitely many commuting flows, as 
indicated in the item (d.l) of Proposition, is assotiated with the arbitrary parameter j in the 
system (140, 141, 146, 147). This remark ends the proof of Proposition. 

5.1 Remark on the eq.(99) 

We will show, that the linear equation (99) is automatically satisfied due to the eq.(119). In 
fact, applying Yla=i G&a* to the eq.(99) and replacing a by a in the result, one gets the next 
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expression for F: 

F\f k) = E E(^ 1 )S n) < ) 4 j;n) - ( 152 ) 

n=l a=l 

The eq.(99) gets the next form: 

U { :; ] = E UWB^KfT^Bffl, n^p (153) 

»,7,7l=l 

in view of the eq.(119), where we have canceled Bp' n ^ assuming that n ^ (3. Eq. (153) is 

Q 

the identity. To show this, we apply •/^'^(i? -1 )^ to the eq.(153), use the eq. (116) to 
simplify the RHS and use the eq.(119) to simplify the LHS of the resulting equation. 

References 

[1] C.S.Gardner, J.M.Green, M.D.Kruskal, R.M.Miura, Phys.Rev.Lett, 19, 1095 (1967) 

[2] V.E.Zakharov and A.B.Shabat, Punct.Anal.Appl., 8, 43 (1974) 

[3] V.E.Zakharov and A.B.Shabat, Funct Anal.Appl., 13, 13 (1979) 

[4] V.E.Zakharov and S.V.Manakov, Punct.Anal.Appl., 19, 11 (1985) 

[5] L.V.Bogdanov and S.V.Manakov, J.Phys.A:Math.Gen., 21, L537 (1988) 

[6] M.J.Ablowitz and H.Segur, Solitons and Inverse Scattering Transform, (SIAM, Philadel- 
phia, 1981) 

[7] V.E.Zakharov, S.V.Manakov, S.P.Novikov and L.P.Pitaevsky, Theory of Solitons. The 
Inverse Problem Method, Plenum Press (1984) 

[8] B.Konopelchenko, Solitons in Multidimensions, World Scientific, Singapore (1993) 

[9] F.Calogero in What is integrability by V.E.Zakharov, Springer, 1 (1990) 

[10] F.Calogero and Ji.Xiaoda, J.Math.Phys, 32, 875 (1991) 

[11] F.Calogero and Ji.Xiaoda, J.Math.Phys, 32, 2703 (1991) 

[12] F.Calogero, J.Math.Phys, 33, 1257 (1992) 

[13] F.Calogero, J.Math.Phys, 34, 3197 (1993) 

[14] F.Calogero and Ji.Xiaoda, J.Math.Phys, 34, 5810 (1993) 

[15] V. E. Zakharov, On the dressing method, in Inverse Methods in Action, edited by P. C. 
Sabatier, Springer- Verlag, Berlin, 602 (1990). 



37 



[16] V. E. Zakharov, Integrable Systems in Multidimensional Spaces, Lecture Notes in Physics, 
153, Springer- Verlag , Berlin, 190 (1982). 

[17] V. E. Zakharov, Multidimensional Integrable Systems in Proceedings of the International 
Congress of Mathematicians, PWN Warsaw, 1225 (1983). 

[18] R. Beals and R. R. Coifman, Scattering, transformations spectrales et equations 
d'evolution nonlineare, I and II at Seminaire Goulaoic-Meyer-Schwartz 1980-1981, exp. 
22, and 1981-1982, exp. 21, Ecole Polytechnique, Palaiseau 

[19] M. J. Ablowitz, D. BarYaacov and A. S. Fokas, Stud.Appl.Math., 69, 135 (1983) 

[20] P.M.Santini, M.J.Ablowitz and A.S.Fokas, J.Math.Phys., 25, 2614 (1984). 

[21] A.Zenchuk, J.Physics A:Math.Gen. 37, nn 25, 6557 (2004). 

[22] A.I.Zenchuk, P.M.Santini, J. Phys. A: Math. Gen. 39 5825 (2006) 

[23] A.I.Zenchuk, math.AP/0603294 

[24] A.I.Zenchuk and P.M.Santini, J.Phys.A:Math.Theor, 40 6147 (2007), nlin.SI/0701031 



38 



